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1 INTRODUCTION 



Gamma-Ray Bursts (GRBs, see e.g. iPiranI |2005|. for a re- 
view) were discovered in the late 1960's ( Klebesadel et all 
1 19731 ). For about three decades, their exact nature re- 
mained elusive. The effort to provide a fast re-pointing of 
the BeppoSAX satellite allowed for the first detection of 
their afterglows at all w avelengths (e.g. I Cost a et al.l 1 19971 : 
Ivan Paradiis et al. Ill997l ). Since then, the specification of 
each generation of instrumentation devoted to GRB studies 
included the need to point to the GRB afterglows quickly. 
Now, with Swift, we have reached the point where we can 
collect data from X-rays to optical within seconds after the 
start of the event, tha nks to fast movi ng telescopes trig- 
gered by GCN notices (jBarthelmvl 1 199^ . This run toward 
fast re-pointing has allowed for prompt optical emission ob- 
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ABSTRACT 

We present the observations of the afterglow of gamma-ray burst GRB 090102. 
Optical data taken by the TAROT, REM, GROND, together with publicly available 
data from Palomar, lAC and NOT telescopes, and X-ray data taken by the XRT 
instrument on board the Swift spacecraft were used. This event features an unusual 
light curve. In X-rays, it presents a constant decrease with no hint of temporal break 
from 0.005 to 6 days after the burst. In the optical, the light curve presents a flattening 
after 1 ks. Before this break, the optical light curve is steeper than that of the X-ray. 
In the optical, no further break is observed up to 10 days after the burst. We failed to 
explain these observations in light of the standard fireball model. Several other models, 
including the cannonball model were investigated. The explanation of the broad band 
data by any model requires some fine tuning when taking into account both optical 
and X-ray bands. 

Key words: gamma-ray: bursts 



serv ations by small auto nomous telescopes for several GRBs 
(e.g. iKlotz et al.ll2009j ). 



Since the mid of the 1990's, 
cal model has emerged in order 



a theoretical canoni- 
to explain the GRB 



phenome non: the fireball model (iRee s fc Mes zaro^ Il992l : 
Mcszaros fc Reesll 19971 : iPanaitescu et a l. 1998). This model 
is based on blast waves emitted by a progenitor during 
black hole formation, that interacts within a jet and/or with 
the surrounding medium to produce all observed emissions. 
These emissions are not thermal, and are proposed to be syn- 
chrotron emission from the electrons accelerated within the 
fireball. This model has the advantages that it is quite easy 
to apply and reproduces, reasonably well, th e observations 
done by BeppoSAX and other telescopes (e.g. iGendre et al.l 
I2OO6I ). The situation has, howev er, radically change d with 
the launch of the Swift satel lite (iGehrels et al.ll2004i ). Swift 
has discovered X-ray fiares jO'Brien et al.l 200d). an early 
fiattening of the X-ray emission (|Nousek et al.n2006l ). and 
seve ral chromatic breaks within the afterglow light curves 
(e.g. iRacusin et al]|2008l ). that were unexpected within the 
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initial standard model. This model has been modified in 
order to explain all these new featu r es, in cluding late en- 
ergy injection ( e.g. [G hiscUini et alj l2007l l. multi-core or 
patchy jets (e.g. IPeng e t al. 2005|, and references therein) , 
and multi-component afterglow (e.g. lWillingale et al.ll2007l ). 
However, these modifications to the standard model were 
made only in order to explain a given feature, and some- 
times lack a multi- wavelength consistency check. In previous 
years, several studies on peculiar events have led the real- 
ization that all available data could not be ex plained due 
to lim itations within the standard model (e.g. iKlotz et al.l 
l2008h . 

In this article, we present observations of GRB 090102 
at optical and X-ray wavelengths, taken over the 10 days 
following the burst. While the optical or X-ray data alone 
could be explained within the framework of the standard 
model, when combining these different datasets, the stan- 
dard model alone cannot explain the observations of this 
burst, indicating that other models need to be tested. The 
paper is organized as follows: we present GRB 090102 in 
Sec. [2] with the data reduction and analysis in Sec. [31 we 
fit several models to these data in Sec. [l] finally, we discuss 
our findings in Sec. [5] before concluding. Errors on spec- 
tral and temporal indices (and the related closure relation 
values) are quoted at the 90% confidence level throughout. 
Reported errors on data points are la errors. 




E 



Figure 1. Image of the field of view of GRB 090102 from 
GROND, taken in the i' band 11 ks after the burst. The GRB 
afterglow is indicated by two ticks. 



2 GRB 090102 



GRB 090102 triggered 
iMangano et~aLl l2009al ) 



Swift-BAT 
on Jan. 2009, 



(trigger 
2""^ at 



338895, 
02:55:45 



UT (hereafter, Tp ) . It w as also recorded by Konus Wind 
(Golenetskii et al' '2009|), and by INTEGRAL/SPI-ACS 
(Marigano ct al. 2009^). The duration of this event, 
as observed by Swift was T90 = 27.0 ± 2.2 s. Its light 
curve presents four p eaks from To — 14 s to To -f 15 s 
ISakamoto et al.ll2009l ). Its time averaged spectrum can be 
represented either using an exponential cu toff power law 
model or a Band model (|Band et al.lll993l ). with parame- 
ters Q = -0.86 ± 0.14, 13 ^ -2.73, and Ep = ^5ltl^s keV 
ixl = 1.075, 64 d.o.f.) and w ith a fiuence of 3.091^-^? x lO"'^ 



erg cm"^ (20 keV - 2 MeV, iGolenetskii et al.ll2009l). Using 



the m easured redshift 
l2009ah . this gives Ei 



(2=1.547, Ide Ugarte Postigo et al.l 
' 5.75 X lO'^^ erg (1 keV - 10 MeV, 
rest frame) and a rest frame Ep of 114911^ keV, fully 
compatible with the so-called Amati relation (|Amati et al.l 
L2OO9, 200i). 

Due to an observational constraint (Earth limb too 
close). Swift slewed to the position of the afterglow with 
a delay, and XRT and UVOT observations began only at 
To -I- 395 s, observing a bright fad ing source in both th e 
XRT and the UVOT fields of view l|Mangano et al.ll2009al ). 
No optical observations were available during the prompt 
phase. The earliest ground search of the optical afterglow 
was performed by the Telescope a Action Rapide pour 
les Objets T ransitoires (TAROT), starting at ro-l-40.8 s 
l|Klotz et al .ll2009a'). The Rapid Eye Mount (REM) robotic 
telescope also obse rved the field starting from ro-^53 s 
(|Covino et al.l |2009| . note however that the start time re- 
ported in the GCN is not correct). A new source was de- 
tected at RA=08''32™58.52", Dec=+33°06'51.1"(J2000.0) 



with an error of 0.3" (see Fig. [T]). This source was then con- 
firmed to be the optical afterglow of GRB 090102. Spectro- 
scopic observa tions were carried by the Nordic Optical Tele- 
scope (NOT) (|de Ugarte Postigo et a"i]|2009ah . the Palomar 
200" Hale telescope ('Cenko et al."2009l), and by the Hobby- 
Eberly Telescope (iCucchiara & Fox 2009,). 

Late afterglow observations were performed by 
the Gamma-Ray burst Optic al/Near-Infrared D etector 
(GROND) 2.50 h after the burst l|Afonso et al.ll2009l ). by the 
Institute de Astroffsica de Can arias 80 cm telescope (IAC80) 
19.2 h after the burst (|de Ugar te Postigo ct al. 2009b), and 
by the P alomar 60-inch tele scope 50 min after the burst and 
later on (|Cenko et al.ll2009l ). The observations continued in 
the follow ing days using the NOT (M alesani et al.ir2009l ) and 
the HST (!l evan et al]|2009h . leading to the detection of the 
host galaxy. 

A radio follow-up was also carried out using the VLA 
and the Westerb ork Synthesis Radio T elescope (WSRT) 
at the 8.46 GHz jChandra fc Fraill l 2009l ) and the 4.9 GHz 
(|van der Horst et al.ll2009l ) wavelengths respectively. No af- 
terglow was detected, leading to upper limits. Finally, the 
MAGIC Cerenkov Telescope also poi nted to the positio n of 
the afterglow, without any detection l|Gaug et al.ll2009l ). 



3 DATA REDUCTION AND ANALYSIS 
3.1 X-ray data 

The data were processed using the ftools version 6.6.2. To 
extract light curves and spectra, a rectangular box and cir- 
cular extraction region were used for window timing and 
photon counting modes respectively, with an extraction re- 
gion of 50 pixels in both cases. XRT observations began 395 
seconds after To in window timing mode. The instrument 
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Figure 2. Light curve of the afterglow of GRB 090102 in X-ray, g', R, i', z', J, H, and K bands from Tables [T] and [2] For clarity, upper 
limits have been omitted. When not plotted, the error bars have the size of the symbols or smaller. These light curves have been corrected 
for the galactic extinction/absorption due to the Milky Way. The X-ray light curve is extracted in the 0.5 to 10.0 keV band. The dotted 
line is the best fit decay law (see text for details). The left axis applies to optical/infrared data, the right axis to X-ray data. 



switched to the photon counting mode fsee lHill et al.|[2005l ) 
at To + 669.1 s, while the count rate was high enough to 
cause pile up. As a matter of consequence, the observation 
between To -|- 669. 1 s and ~ To + 3000 s suffers from severe to 
moderate pile up. This area was discarded in the following 
spectral analysis, a nd corrected within th e light curve using 
the prescriptions of IVaughan et all (|2006l '). 

All spectra were rebinned in order to include 20 net 
counts within each bin, and fit using the statistic. The 
spectral model was composed of a power law continuum 
absorbed by our own g alaxy {Nh = 4.1 x 10^° cm~^, 
iDickev fc LockmanI Il990l ) and by the host galaxy at z = 
1.547. Extracting spectra at several epoch, we did not find 
any evidence for spectral variation during the whole follow- 
up of the afterglow. The spectrum is well modeled {xt = 
0.92, 102 d.o.f.) by a single power law with spectral index 



Px ~ 0.83 ± 0.09 absorbed in our own galaxy and by extra- 
galactic absorbers {Nh = S.Slj'j x 10^^ cm~^ when located 
within the host galaxy). 

The hght curve was extracted within the 0.5-10.0 keV 
band and rebinned in order to obtain at least 25 counts 
per bin. All decay index indicated below are derived from 
fits using the statistic. We did not observe strong flares 
(see Fig. A single moderate flare can be seen in the un- 
binned hght curve during orbit 2 (corresponding to the pe- 
riod ~ 2000 — 4000s after the trigger), which has been ex- 
cluded from the temporal analysis. The complete hght curve 
can be adequately flt {xt = 1.18, 65 d.o.f.) using a single 
power law with a decay index ax = 1-34 ± 0.02. A bro- 
ken power law can also represent the data (xt ~ 0.99, 63 
d.o.f.) with Qi.x = 1.29 ± 0.03, 02. x = 1.48 ± 0.10, and 
a break time tb = ISTOO^ggpg" s. With the observed steep- 
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ening, Aa = 0.19 ± 0.11, this break could be interpreted 
as the cooling frequency passing through the observation 
band. However, this should be associated with a spectral 
break not supported by the data. An F-test check on the 
break existence gives a probability of 0.26, i.e. a value not 
conclusive. One may note however that in this special case 
(large error on tb that make it compatible with within 
3(t), the probability deri ved from the F- test should not be 
considered as valid (see iProtassov et al.l [2002i ) . As a mat- 
ter of consequence, in the following we will consider both 
hypotheses (single power law or broken power law), and re- 
port early and late X-ray data as data taken before and after 
the temporal break respectively. 



3.2 Optical data 

At the position of the afterglow, the G alactic extinc- 
tion i s E{B - V) ^ 0.047 according to ISchlegel et all 
We corrected all magnitudes using this estima- 
tion. Assuming Rv=3.1, this gives in particular ylR=0.14 
|Peilll992l ). The specific data processing relative to all in- 
struments is indicated in the following subsections. We 
cor npleted this sampl e by using reported observations 



bv ICenko et al.l J 2OO9I): 'd e Ugarte Postigo et al.] (|2009bh : 



iMalesani et al.l ((2009). This extended sample of data taken 
in the R band is listed in Table [T] 



3.2.1 TA ROT data 



We used data from TAROT Calern (|Klotz et al.ll2009dr ) that 
started an exposure of the field of GRB 090102 at To-|-40.8 s 
(duration 60 s) with the tracking speed adapted to obtain a 
small ten pixel trail. This technique was used in order to ob- 
tain continuous temp oral information during the exposure 
fe.g. lKlotz et al.l|2006l ). The spatial sampling was 3.29"pix~^ 
and the FWHM of stars (in the perpendicular direction of 
the trail) was 2.05 pixels. Only the first exposure was per- 
formed with this technique. Successive images were tracked 
on the diurnal motion using exposure times increasing from 
30 s to 180 s. Images are not filtered. 

We used the star USNO-Bl 1231-0206902 as a constant 
template for calibrat ion during the tr ailed exposure. Con- 
trary to the claim of iKlotz et al] (|2009a '). TAROT did not 
observe a rise during the first part of the trailed image. In- 
deed, the profile is convolved by the Point Spread Function 
(PSF) that induces a slope at the sta rt and at th e end o f 
the trail. The rapid analysis done by iKlotz et al.l (|2009al ) 
confused this slope with a rise at the start of the GRB af- 
terglow trail. This has been corrected for by using the flux 
of the template (constant) star. 

In order to precisely cross-calibrate the observations 
performed with REM and TAROT, we used several USNO 
stars as relative standards. All flux measurements were done 
using standard aperture photometry with the AuDela soft- 
war<|3. USNO-Bl 1231-0206902 has been used to build the 
PSF in TAROT images in order to fit the flux of the GRB. 
TAROT observations are not filtered. We thus converted 
all TAROT instrumental magnitudes to standard R magni- 
tudes, using the following formula: R= C — 2.5 x log{F) + 

^ http://www.audela.org/ 



(V — R) X Qfr. In this equation, ar — 0.4 is determined using 
calibration stars with various (V-R) values located within an 
open cluster. C is calculated for each image from the PSF 
flux, the R and (V-R) measurements of the reference stars. 

3.2.2 REM data 

Early time optical and NIR data were collected using the 60 
cm robotic telescope REM (|Zerbi et al.ll200ll : ICovino et al.l 
[2004) located at the ESO La Silla observatory (Chile). The 
telescope simultaneously feeds, by means of a dichr oic, the 
two focal instruments: REMIR l|Conconi et al.ll2004l ') a NIR 
camera, operating in the range 1.0-2.3 fim (z', J , H and K), 
and ROSS (REM Optical Slitless Spectrograph. iTosti et"al] 
l2004h an optical imager with spectroscopic (slitless) and 
photometric capabilities (V, R, I). Both cameras have a field 
of view of 10 X 10 '. 

REM reacted promptly and began observing 103 s after 
the GRB trigger time, following the event for ~1 hour (see 
Table [T]). For the first 500 s REMIR observations were per- 
formed only in the H filter with increasing exposure times, 
then all the NIR filters were used in sequence. In the opti- 
cal range only R band images were taken, for a total of 38 
consecutive images. The transient was only detected up to 
500 s both in the R and H bands, after which its luminosity 
falls below the instrument detection limits for all filters. 

Each single NIR observation with REMIR was per- 
formed with a dithering sequence of five images shifted by a 
few arcseconds. These images are a utomatic ally elaborated 
using the proprietary routine AQuA (|Testa et al.:..2004 '). The 
script aligns the images and co-adds all the frames to obtain 
one average image for each sequence. 

Each R image was reduced using standard proce- 
dures. A combination of the ESO-MIDaI I DAOPHOT task 
(Stetson, 1988 ) and Sextractoi[f| package ( Bertin &: Arnoiiti 
Il996h was used to perform standard aperture photometry. 
The photometric calibration for the NIR was accomplished 
by applying average magnitude shifts computed using the 
2MAS^| catalogue. 

The optical data were calibrated using instrumental 
zero points and checked wit h observations of standard stars 
in the SA92 Landolt field (|Landolt|[l993 'l. Results are dis- 
played in Tables [T] and [21 



3.2.3 GROND data 

The Gam ma-Ray Burst Opt ical/Nearlnfrared Detector 
(GROND, iGreiner et al] l2008t ) mounted at the 2.2 m 
ESO/MPI telescope at La Silla observatory imaged the 
field of GRB 090102 simultaneously in g'r'i'z' JHK starting 
9.01 ks after To under clear sky conditions ((Afonso et al.l 
2009). At that point the location of the burst was becoming 
visible above the pointing horizon of the telescope. A total 
of 12 g'r'i'z' images with integration times of 66 s, 115 s and 
375 s were obtained during the first night post burst. At the 
same time, NIR frames were acquired with a constant expo- 
sure of 10 s. In addition, the field was observed by GROND 



^ http:/ /www. eso.org/projects/esomidas/ 

^ http:/ /astroa. physics. metu.edu.tr/MANUALS/sextractor/ 

* http: / / www.ipac.caltcch.edu / 2mass / 
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Table 1. Table of data taken in the R band, corrected for galactic extinction. GROND data in r' band are transforme d to Vega R 
system . Errors are given at Itr level. Data from Palomar, IAC80 and NOT are taken from iCenko et al.i ()2009i" l : idc Ugartc Postigo et al.l 
ll2009d) : [Malesani et al] | |2009^ respectively. 



Mid time 


Exposure 
duration 


Magnitude 


Instrument 


Mid time 


Exposure 
duration 


Magnitude 


Instrument 


(s) 


(s) 






(s) 


(s) 








o.u 


1 n -1- fi 9 
lo.u m u.z 


TA ROT 
1 l\lX\J 1 


yzoD 


DO 


90 A'^ -\- n 0'^ 
ZU.40 m u.uo 






u.u 


J-O.O ^ U.Z 


TAROT 


yo ( o 


uo 


90 AS -k Ofi 
ZU.4:0 ^ u.uo 




uO.cS 


D.U 


1 Q Q -U fl 9 
lo.o it U.Z 


TA ROT 


Qc;i A 
yOi4 


1 1 

110 


on AA -L r\ rxA 

ZU. 44 It U.U4 




OJ-.o 


o.u 


lO.'rt III U.Z 


TA ROT 


Q71 1 
y J 11 


110 


on cri -1- 09 
ZU.Ol III U.UZ 




f\'7 S 


o.u 


1 ft -1- n 9 

lo.o III u.z 


TA ROT 


yyuo 


110 


on CO _|_ n 09 

zu.DZ zn U.UZ 




( O.O 


O.U 


lo.o it U.z 


TA ROT 


QOVQ O 

yy / o.y 


1 90 
IZU 


on AA -\- r\ 07 
ZU.44 it U.U/ 


Palomar 


/ y.cs 


O.U 


1 /I n -1- n 9 

14. U it U.Z 


TA ROT 


1 0001 

luuy 1 


1 1 
110 


on c;n -U n nQ 
zU.oU it U.Uo 




oo.o 


ou 


14.10 It u.iu 




1U4ZO 


0/0 


on cc _|_ n 09 
ZU.OD It U.UZ 




fSo.o 


O.U 


1/11 _!_ n 9 

14.1 It U.Z 


TA ROT 


lUoDo 


O / 


on fi'i -I- n ni 

ZU.OvS it U.Ul 




Q1 54 
yi.o 


o.u 


1 A -1- n 9 

14. It U.Z 


TA ROT 


1 1 ^90 


0/0 


on fif^ -)- 01 
ZU.DO It U.Ul 




Q7 S 


o.u 


14.4 It U.O 


TA ROT 


I 1 7S0 

II ( oU 


0/0 


on 71-1-0 09 
ZU. / 1 III U.UZ 




122 7 


29 4 


1 A 7 -1- n 1 

I ^ U. 1 


TAROT 


11911 5 


120 


OCX f,p, J- n r\Q 
ZU.OU ^ U.uo 


Palomar 




OU 


1 A QQ _|_ n 90 

14. yo It u.zu 


rxWoo 


14UOO.O 


1 90 
IZU 


on QQ -U OQ 

zu.oy It u.uy 


Palomar 




OU.O 


1 9 -1- 9 
10. z ^ u.z 


TA ROT 


10144. i 


1 90 
IZU 


on Q7 -1- OQ 

zu.y 1 It u.uy 


Palomar 




oU 


1 c; 1 Q _i_ n oo 
10. lo it U.zU 


ROQQ 


1 fi 1^99 R 
loOZZ.O 


DUU 


91 11 -U n Of^ 
zi.ii it u.uo 


Palomar 


lyD.o 


oU.O 


1 c; T' _!_ n o 
10. ( it U.z 


TA ROT 


zUDoi.O 


DUU 


oi Qn -U n n(^ 
zi.oU it u.uo 


Palomar 


909 9 


OU 


1 c; CO _|_ fi on 
10. Oo It u.zu 


rtW oo 


ZZ / OD. 1 


fiOO 
DUU 


91 Afi -1- 07 
Z 1 .4D It u.u / 


Palomar 


241.1 


30 


15.78 lb 0.20 


ROSS 


24832.6 


600 


21.74 ± 0.08 


Palomar 


252.3 


67.8 


16.0 lb 0.2 


TAROT 


27994.7 


1200 


21.60 ± 0.05 


Palomar 


280.8 


30 


15.98 lb 0.30 


ROSS 


32188.5 


1200 


21.90 ± 0.06 


Palomar 


358.6 


90 


16.48 lb 0.30 


ROSS 


36411.3 


1200 


21.81 ± 0.11 


Palomar 


402.3 


185.4 


16.8 lb 0.2 


TAROT 


69120.0 


1800 


22.48 ± 0.30 


IAC80 


619.5 


330 


17.48 lb 0.40 


ROSS 


81566 




22.98 ± 0.20 


NOT 


943.8 


685.2 


17.9 lb 0.3 


TAROT 


97805 


1500 


23.16 ± 0.07 


GROND 


2283.0 


2173.2 


19.1 lb 0.5 


TAROT 


99200 




23.38 ± 0.30 


1AC80 


2328.5 


1380 


19.38 lb 0.60 


ROSS 


99600 


1500 


23.23 ± 0.10 


GROND 


3190.2 


120 


19.63 lb 0.13 


Palomar 


182130 


1500 


23.55 ± 0.12 


GROND 


5266.2 


120 


19.90 lb 0.06 


Palomar 


183448 




23.58 ± 0.20 


NOT 


6014.4 


5275.2 


19.9 lb 0.4 


TAROT 


183945 


1500 


23.67 ± 0.11 


GROND 


7704.3 


120 


20.29 lb 0.07 


Palomar 


185758 


1500 


23.67 ± 0.11 


GROND 


9042 


66 


20.42 lb 0.04 


GROND 


187575 


1500 


23.57 ± 0.09 


GROND 


9158 


66 


20.45 lb 0.03 


GROND 


264553 




23.92 ± 0.09 


NOT 



at days 1 and 2 after the trigger. D ata reduct ion was per- 
formed using standard IR AF tasks jTodv"l993), similar to 
the procedure outlined in iKruehler e t al. (2008). Absolute 
photometry of the afterglow was measured relative to field 
stars from the SDSS and 2MASS catalogs, which yields av- 
erage accuracies of 0.04 mag in g'r'i'z' and 0.07 mag in JHK, 
and are reported in Tables [T] and [21 

3.2.4 Temporal and spectral behaviors 

The optical light curve can be fit using a broken power law 
model, with a break time on the order of 1000 s. As can be 
seen in Fig.[2l the hght curve flattens at that time. It is worth 
noting that this flatt ening cannot be du e to the host galaxy, 
fainter than 24 ijLevan et al ] |2009l ) when the afterglow 
magnitude is of the order of R~ 19 — 20. Before the break, 
the decay index is ai,o ~ 1-50 ± 0.06, steeper than the 
X-ray decay index at the same date (see Table After the 
break, the decay index is 02,0 = 0.97 ±0.03, flatter than the 
X-ray decay index at the same time. We do not observe any 
color variation (both in optical and in X-rays) before and 
after the temporal break seen in optical (see Table |4|. If we 
assume the host extinction to be null, the value of R — H = 



Table 3. Summary of the properties of GRB 090102. 



Property 


GRB 


090102 


Early optical decay index 


1.50 


± 


0.06 


Late optical decay index 


0.97 


± 


0.03 


X-ray decay index 


1.34 


± 


0.02 


X-ray spectral index 


0.83 


± 


0.09 



2.4 implies Po,obs = 1.32. A non zero extinction would cause 
a flatter Po, thus one should consider this as an upper limit: 
Po < 1.32. The value of the optical-to- X-ray spectral index 
is Pox > 0.53 ± 0.04. 

3.2.5 Optical extinction 

The estimation of the optical extinction is a non-trivial prob- 
lem, as we need to assume a spectral slope of the after- 
glow continuum in order to deduce the amount of extinction 
needed in the optic al to reproduce that continuum (see e.g. 
IStratta et al-llioO^ l. In this article we prefer to focus on the 
constraints we can set on the GRB models, and thus we can- 
not make any assumption on the afterglow continuum. As a 
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Table 2. Table of optical/infrared data taken in other bands. GROND data in g', i', z' bands are transformed to Vega system. All 
magnitudes are corrected for local extinction due to the Milky Way, and errors are given at la. 
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Table 4. GRB 090102 afterglow R-H color. All magnitudes are 
corrected for local extinction due to the Milky Way. 



Mid time 


R-H 


(s) 


color 


88.2 


2.42 ± 0.25 


102.9 


2.47 ± 0.33 


152.1 


2.53 ± 0.24 


166.8 


2.25 it 0.23 


197.1 


2.48 ± 0.27 


11124 


2.39 ± 0.08 



Table 5. GRB 090102 optical SED fit result. We indicate for 
various models the best fit parameters and the goodness of the 
fit. Within the fit, the jSp parameter was not allowed to become 
negative. In all cases, the redshift was fixed to 1.547. 



Model 



\ host 



/3o 



xi 



Milky Way 
Large Magelanic Cloud 
Small Magelanic Cloud 

Calzctti law 



0.35 ± 0.29 

"S2+''-23 
"•°^-0.3S 

0.04^ 



0.77 



+0.13 
0.04 
0.23 
0.77 



52+°'^2 
'^■''^-0.36 

og+^-so 

^•^''-0.28 

003+013 

U.»0_g 22 

1+1.0 
-0.1 



o.r 



0.07 
0.07 
1.35 
1.19 



d.o.f. 



matter of consequence, we have tried to fit the optical-near 
infrared SED without using the X-ray information, assuming 
a simple power-law modej^. We have tested the most com- 
mon models used, i.e. the Milky Way, the Large and Small 

® One should think that this is already an hypothesis on the data, 
i.e. that there is no spectral break located between the g' and the 
K bands. 



Mage l lanic Clouds and th e Calzetti law (see ICalzetti et al] 
I1994I : IStratta et al.l |2004| . and references therein) . Results 
are reported in Table (5] As can be seen from that Table, all 
models are acceptable. However, the small value of xi re- 
ported in some case indicates a statistical problem. In fact, 
the spectral index and the host extinction are degenerate 
parameters, as can be seen from Fig. 3] if the errors are de- 
fined in Table [5] this is only due to the fact that both A„ 
(for physical reason) and /3o (by hypothesis) are not allowed 
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Frequency (Hz) 



Figure 3. Spectral energy distribution of the afterglow of GRB 
090102 observed 11.3 ks after the burst. Frequencies are expressed 
in the observer frame. The solid line represent the best fit X- 
ray model. Dotted line is the direct extrapolation of this model 
toward optical frequencies (assuming no extragalactic extinction). 
As there may be an X-ray flare at that time, we also rescaled this 
extrapolation (dashed line) to the level of the X-ray continuum 
at the same time. 



and cannot be used to derive an estimation of the continuum 
flux density. Assuming that the optical and X-ray emissions 
are produced by the same mechanism, the true value is lo- 
cated between the data points and the dotted line on Fig. 
[3l corresponding to an extinction range of - ~ 2.4 mag. 



4 MODEL FITTING 




4.1 Standard model 



0.6 0.8 
Spectral index | 



We fi rst used the prescriptions of IPanaitescu fc Kumar! 
(|200d ) and the X-ray data alone in order to test th e standard 
model. They are summarized in Tab le [6] (from ISari et al] 
1 19981 : ISari fc Piranll 19991 : lRhoads|[l997l ) . As one can see from 
this Table, in order to accommodate the standard model 
and the measurements made in X-ray, we should have the 
cooling frequency located above the X-ray band during 
the entirety of the observation. In such a case, a value of 
p = 2/3x = 2.7 ± 0.2 can fit both the spectral and tempo- 
ral X-ray indices if the surrounding medium has a constant 
density. However, a constant density medium implies that 
the cooling frequency is decreasing as f""'^. Any crossing 
of this frequency within the X-ray band would appear as a 
simultaneous increase in the spectral index by 0.5 and in 
the temporal decay index by 0.25. While the latter could 
be accommodated by the broken power law model, the for- 
mer is not supported by the data. This implies that the 
cooling frequency should be above the X-ray band up to 
~ 5 X 10^ s after the event, or above 800 keV at the start 
of the X-ray observation (thus implying very small values of 
tB and/or Stot, the energy fraction located within the mag- 
netic field and the total energy of the fireball respectively 
IPanaitescu fc Kumar 2000). This is unlikely; moreover, in 
such a case, according to the model, the optical light curve 
should decay either with the same slope than the X-ray light 
curve or with a decay index of 0.5 (Panaitcscu fc Kum_3 
I2OOOI) : both hypotheses are ruled out by the data. 

In conclusion, the standard model alone, in its simplest 
J 2 expression, cannot reproduce the global spectral and tem- 
poral behavior of the afterglow of GRB 090102. 



Figure 4. Contour plot in the fio-Ay plane for a Large Mag- 
ellanic Cloud extinction model. Contours are indicated for the 
68%, 90% and 99% confidence levels. As one can see from the 
pronounced "banana-shape" of the plot, there is a strong degen- 
eracy between these two parameters. 

to become negative. As a matter of consequence, we cannot 
derive the host extinction from the optical photometric data 
alone. 

However, we are able to provide two comments. First, 
in the host frame, we have an X-ray measurement of Nh = 
8.812 5 ^ lO^^ cm~^. This may imply that significant local 
extinction is present. Within our galaxy, this would lead to 
Ay = 4.9l}'4, th at can be c onsidered as an upper limit (see 
the discussion of IStratta et a l. 2004). 

Finally, in the host frame we have measurements rang- 
ing from 180.1 nm to 852.2 nm, i.e. covering the UV to ~ I 
bands. This last band should be the least affected by the host 
extinction. This is, however, not a band free of extinction. 



4.2 Two component models 

The afterglow of GRB 061126 features the same tem- 

J ioral and spectra l behavior as the on e of GRB 090102 
Perlev et al] I2OO8I : [Goinboc et all |2008| . see also Table [3]). 
These authors used a multi-component model rather than 
the standard fireball model in order to explain the broad- 
band data of this burst. We thus checked if inserting an 
additional component to the fireball model could reproduce 
the data of the afterglow of GRB 090102. 

The first step was to decide in which band we see the sig- 
nature of the extra-component, and in which band the fire- 
ball model alone can reproduce both the spectral and tem- 
poral properties. In the standard fireball model, the fact that 
the optical and X-ray light curves do not present the same 
decay implies that a specific frequency (either t he cooling 
or th e injection frequency, see their definition in iMeszaroa 
I2OO6I ) is located between the optical and X-ray bands. Be- 
low the injection frequency, the decay index should be 0.5. 
With an optical decay of Ofi.o = 1-50 ± 0.06, we can rule 
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Table 6. Closure relationships in the standard fireball model computed using the spectral and temporal information from the X-ray 
band. The specific frequency is i/c (the cooling frequency) in the slow cooling regime, or Um (the injection frequency) in the case of fast 
cooling. All errors are quoted at the 90 % confidence level. Due to the absence of measured spectral index in the optical, optical values 
of the relationships are not reported 



Medium class 


Cooling 


Specific 




Closure 






Observed value 




and 


regime 


frequency 




relationship 




X-rays (0.5-10.0 keV) 




geometry 




position 








(single power 
law hypothesis) 


(early broken 
power law) 


(late broken 
power law) 


Isotropic Wind 


Fast 


l^m < l^X 


a 


- 1.5/3 + 0.5 = 





0.6 + 0.2 


0.5 + 0.2 


0.7 + 0.3 






Vm > VX 




a - 0.5/3 = 




0.9 + 0.1 


0.9 + 0.1 


1.1 + 0.2 




Slow 


Vc < I^X 


a ■ 


- 1.5/3 + 0.5 = 





0.6 + 0.2 


0.5 + 0.2 


0.7 + 0.3 






Vc > Ux 


a 


- 1.5/3 - 0.5 = 





-0.4 + 0.2 


-0.5 + 0.2 


-0.3 + 0.3 


Isotropic ISM 


Fast 


I'm < l^X 


a 


- 1.5/3 + 0.5 = 





0.6 + 0.2 


0.5 + 0.2 


0.7 + 0.3 






Vm > l^X 




a - 0.5/3 = 




0.9 + 0.1 


0.9 + 0.1 


1.1 + 0.2 




Slow 


I^c < Vx 


a 


- 1.5/3 + 0.5 = 





0.6 + 0.2 


0.5 + 0.2 


0.7 + 0.3 






> t^X 




a - 1.5/3 = 




0.1 + 0.2 


0.0 + 0.2 


0.2 + 0.3 


Jetted fireball 


Slow 


Vc < Vx 




a - 2/3 = 




-0.3 + 0.2 


-0.4 + 0.3 


-0.2 + 0.3 






Vc > Vx 


a 


- 2/3 - 1.0 = 





-1.3 + 0.2 


-1.4 + 0.3 


-1.2 + 0.3 



out this hypothesis, and deduce that the cooling frequency 
needs to be located between the optical and X-ray bands. 
This is not compatible with the X-ray data (see Table [6]): 
an extra component in the X-ray band is the only solution 
to solve this issue. 



Under the above assumption, the optical band can be 
described by the fireball model and it is the X-ray band that 
needs an extra component in addition to the fireball in order 
to be described. This hypothesis needs to be confirmed by 
checking if the optical band can indeed be fit with the stan- 
dard model. In fact, the fiattening of the optical light curve 
around 1 ks is not possible in the standard model. However, 
several "common additions" to the standard model can pro- 
vide an explanation: the addition of a reverse shock in the 
early optical data or a change in the surr ounding medium 
density profile (the termination shock, see IChevalier et al. I 
l2004l : Ir amirez- Ruiz et al. II2001I ). In the latter case, we can 
derive a value of p = 4/3a2,o + 1 = 2.29±0.04. The physical 
position of the termination shock cannot be constrained due 
to the presence of a second component in the X-ray data. We 
can thus conclude that our hypothesis is valid: the optical 
band can be fit with the standard model. 



Lastly, one then needs to provide for the physical na- 
ture of the additional X-ray component. The requirements 
are that it needs to be: non-thermal, at all times brighter 
than the forward shock emission responsible for the optical 
emission, decaying monotonously without spectral variation, 
and producing few (if any) optical emissions. We do not con- 
clude on the mechanism that produces these extra photons 
in X-ray. However, the latter requi rement is the more con- 
straining one, as already noted by IPerlev et al.l l|2008h for 
GRB 061126: it is quite difficult to accommodate a mech- 
anism that produces X-rays without also producing optical 
photons. Note additionally that late in the light curve (after 
1 ks), the standard model implies that the X-ray spectral 
index is p/2 ^ 1.2. The additional component needs to pro- 
duce hard photons in order to obtain Px ~ 0.83 ± 0.09. 



4.3 The cannonball model 

In the cannonball model, GRBs and their afterglows are pro- 
duced by bipolar jets of highly relativistic plasmoids (CB) of 
ordinary matt er ejected in the b irth of neutron stars or black 
holes (see e.g. iDado et al.|[2009 l. and references therein). 

A prompt, hard X-ray /gamma-ray pulse is produced by 
the thermal electrons in the CB plasma, via inverse Comp- 
ton scattering of photons emitted/scattered into a cavity 
created by the wind/ejecta blown from the progenitor star 
or a companion star long before the GRB itself. Slightly 
later, when the CB encounters the wind/ejecta, the elec- 
trons of the ionized gas in front of it are swept into the CB 
and are Fermi accelerated by its turbulent magnetic field, 
emitting synchrotron radiation that dominates the prompt 
optical emission and the broad band afterglow emission. 
Thus, the observed emissions are due to two components: a 
short lasting Inverse Compton Scattering component, seen 
only at high energy, and a classical synchrotron emission 
( Dado fc Par.. 2009 ). 

F ollowing the notations and prescriptions of lDado et al.l 
(|2009l ), we can explain all the X-ray data assuming that we 
observe the synchrotron emission of a CB traveling within 
a medium of constant density (the so-called ISM) with the 
bend frequency located below the X-ray band. From the 
model, which gives, in that case fix ~ p/2, we can derive 
p = 1.7 ± 0.2. This value of p implies ax = P + 0-5 = 
p/2 + 0.5 = 1.35 ±0.1, in perfect agreement with our finding 
[ax = 1-34 ± 0.02). We note that this value of p is lower 
than 2, and thus there should be a break/cutoff located at 
high energy so that the electron energy remains finite. 

In the optical , things are more complicated. Accord- 
ing to iDado et al.l (|2009), the spectral index in the opti- 
cal should evolve from /3o ~ 0.5 during the prompt emis- 
sion to Po ~ 0.83 ± 0.09 (i.e. the same as that in X-ray) 
during the late-time afterglow. Both values are compatible 
with the constraint derived from the color of the afterglow 
(P ^ 1.32). During the prompt optical emission, the de- 
cay index should be ao = /3o + 1 ~ 1-5 (see Eq. 33 in 
iDado et al.ll2009l ). This is consistent with the observed value 
of ai,o = 1-50 ± 0.06. The flattening of the optical light 
curve may be explained by a transition from a wind envi- 
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ronment to an ISM (as in the fireball model, see above). In 
such a case, the light curve should reach an asymptotic de- 
cay with Q2.0 = ax = 1.34 ± 0.09 (see Eqs. 29 and 30 in 
iDado et al.ll2009l ). We observed a value of 02,0 = 0.97±0.03. 
However, the value of 02,0 = 1-34 ± 0.09 is only an asymp- 
totic value, and the observed could vary if the observation 
is close to the time of the crossing of the bend frequency. 
In fact, we should observe an effective temporal decay index 
a2{t) = l3oit) + l/2 evolving from 02 ~ 1 to ^2 = 1.34±0.02. 
Being close to the asymptotic value, the measurement is 
marginally compatible with the CB model expectation. Lack 
of late time temporal and spectral data prevents a confirma- 
tion of this prediction. Indeed, once reaching the asymptotic 
values, the CB model indicates that we should observe the 
same spectral and temporal behavior for the optical and X- 
ray bands. Looking to Fig. [S] this is hardly supported by 
the data. However, a large gray extinction, with A/ ~ 2.5 
would make compatible the optical and X-ray emissions. In 
the absence of information about the host extinction, we 
cannot rule out such an extinction, and the CB model may 
indeed reproduce all the data with the expense of a few fine 
tunings. 



5 DISCUSSION AND CONCLUSIONS 

We have presented observations of the afterglow of 
GRB 090102. We used data taken by optical telescopes and 
by the XRT on board the Swift spacecraft. The light curve 
of this event features an unusual behavior. In X-rays, it 
presents a very smooth light curve with no hint of a tempo- 
ral break. In the optical, the light curve presents a steep-flat 
behavior, with a break time at ~ 1 ks after the burst. This 
unusual light curve was also observed with GRB 061126. The 
standard fireball model cannot reproduce the observations 
without the addition of several components. The addition of 
a component at high energy solves the discrepancies between 
the optical and X-ray emissions. The optical data could then 
be interpreted as either due to a termination shock, locating 
the end of the free- wind bubble at the position of the optical 
break; or as a normal fireball expanding in an ISM, with a 
reverse shock present at an early time (less than 1 ks af- 
ter the event). Any radio observations that could constrain 
the ISM density value would also constrain the wind den- 
sity, and thus a key property of the progenitor of this GRB. 
However, the weak point of this model is that requires strong 
fine tuning in order to suppress the optical emission of the 
additional component, responsible for the X-ray emission. 
The cannonball model can also partly reproduce the data. 
It appears clear, however, that in order to explain the broad- 
band emission, some fine-tuning of this model is mandatory, 
likewise for the fireball model. 

The strange behavior of GRB 090102 provides a useful 
test-event for models that attempt to reproduce and explain 
the GRB phenomenon. The choice of the models tested here 
was strongly limited by the very small number of models 
that have published formulae of the temporal and spectral 
indices for GRB afterglows. Of course, we do not pretend to 
be exhaustive in our testing as the array of alternative mod- 
els is huge. Swift has shown the diversity of GRB afterglows, 
and it may now be time for observers to test and reduce the 
diversity of GRB afterglow models. 
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